Failure of standard approximations of the exchange couphng in nanostructures 
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We calculate the exchange coupling for a double dot system using a numerically exact technique 
based on finite-element methods and an expansion in 2D Gaussians. Specifically, we evaluate the 
exchange coupling both for a quasi-one and a two-dimensional system, also including an applied 
magnetic field. Our numerical results provide a stringent test of standard approximation schemes 
(e.g., Heitler-London, Hund-Mulliken, Hubbard), and they show that the standard methods do not 
have reliable predictive power even for simple model systems. Their value in modeling more realistic 
quantum-dot structures is thus cast in serious doubt. 
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I. INTRODUCTION 

The possibility of coherent manipulation of electron 
spins in low-dimensional nanostructures, aimed at future 
large-scale quantum information processing,— calls for a 
thorough understanding of the spin interactions at play. 
In the proposal for quantum computing with quantum 
dots by Loss and DiVincenzo, the exchange coupling be- 
tween the spins of electrons in tunnel-coupled quantum 
dots was envisioned as the controllable mechanism for co- 
herent manipulation of spin qubitsiii^ Recently, this fun- 
damental building block of a possible future solid-state 
quantum computing architecture was realized in an ex- 
periment, demonstrating electrostatic control of the ex- 
change coupling.— 

In this paper we present numerically exact finite- 
element methods for calculations of the exchange cou- 
pling between electron spins in tunnel- coupled quasi- 
one and two-dimensional quantum dots. Such struc- 
tures have already been under intensive theoretical 
investigation using various numerical methods, e.g. 
based on an exact diagonalization of the underlying 
Hamiltonianii^i^ii^i^ii^ or using quantum-chemical ap- 
proaches like self-consistent Hartree-Fock methods.— 
Such numerical approaches often require extensive nu- 
merical work. Therefore, much attention has been de- 
voted to simple approximations which lead to closed-form 
analytic expressions for the exchange coupling.— '^^d^ It is, 
however, not immediately obvious to what extent these 
approximations yield correct predictions, and where they 
break down. For example, in a recent worki^ the validity 
criterion for such approximations was the requirement 
that the exchange coupling at zero magnetic field always 
must be positive. A criterion like this can only provide 
a necessary condition for an approximate scheme to be 
acceptable. 

The aim of this work is to provide a quantitative com- 
parison of the Heitler-London, the Hund-Mulliken, and 
the Hubbard approximations, applied to a simple model 
potential of a double quantum dot, with numerically ex- 
act results. In particular, we focus on the case, where 



the distance between the two quantum dots is short, 
such that the single-dot electron wave functions have a 
large overlap. For short distances, the exchange cou- 
pling can reach values on the order of several meV, mak- 
ing it sufficiently large to exploit and observe in exper- 
iments, and our comparative study is thus highly rele- 
vant for on-going experimental activities within the field. 
The finite-element methods used here allow an easy im- 
plementation using available numerical packages;^ also 
when finite magnetic fields are included, which strongly 
infiuence the exchange coupling in two-dimensional ge- 
ometries. We find that the approximative schemes may 
provide reasonable predictions of the exchange coupling 
for certain parameter ranges, while they fail, also qual- 
itatively, for short distances, even for the simple model 
potential considered here. Their value in modeling more 
realistic quantum-dot structures used in experiments is 
thereby cast in serious doubt. 




FIG. 1: (Color online). Double quantum dot and numeri- 
cally calculated charge density. The double quantum dot is 
described by the potential V(v) given in Eq. (|4]) (here with 
Of = 1, hujo = 4 meV and d/ro = 1). The two-dimensional 
contour plot shows how the charge of two electrons in a sin- 
glet spin state is distributed within the double quantum dot. 
With finite tunnel-coupling between the two quantum dots, 
the spins of the electrons interact due to the exchange inter- 
action. 
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II. DOUBLE QUANTUM DOT MODEL 

Experimentally, electrons can be confined in double 
quantum dots using metallic gates on top of a semi- 
conductor heterostructure^ii^ii^ or across a nanowireiii^ 
or a nanotubej ^^i^° By suitable electrostatic gating such 
Coulomb-blockade double quantum dots can be brought 
into a few-electron regime;^ where only a single electron 
occupies each of the two quantum dots. In this regime, 
the spin and charge dynamics are described by a two- 
electron Hamiltonian of the form 



where 



i7(ri, rs) = h{ri) + hir^) + C(|ri - raj). 



C{\r,-r,\) = 



47rereo ri - r2 



(1) 



(2) 



is the Coulomb interaction and the single-particle Hamil- 
tonians are 



spin Hamiltonian, 7i = J 81-821^ The task is to calculate 
the exchange coupling J as function of various parame- 
ters, e.g., the distance between the quantum dots and 
the applied magnetic field. A magnetic field only affects 
the exchange coupling significantly in two-dimensional 
geometries and we consequently concentrate on the in- 
clusion of a magnetic field in the two-dimensional case 
a=l. 



III. VALIDITY OF APPROXIMATE METHODS 

A. Quasi-one dimensional limit 

We first consider the quasi-one dimensional limit a ^ 
1, which may be relevant, e.g., for describing confined 
electrons in nanowires. In this limit we integrate out the 
motion in the y-direction and consider an effective one- 
dimcnsional model reading 



(3) 



H = h{xi) + h{x2) + Ca{\xi - X2\), 

where the single-electron Hamiltonian is 



with V{r) denoting the confining potential. As in many 
realizations of double quantum dots wc assume that the 
motion of the electrons is restricted to maximally two 
dimensions, i.e., r = {x,y). The inclusion of a magnetic 
field is discussed below. 

As an illustrative example^^ we consider a simple dou- 
ble dot potential reading^i^ 



h{x) 



2m 



Vix), 



nx) = ^[min{(x-d)^(.T + d)n] 



and we have introduced 
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V{r) = ^ [min{(x - d)^ + d)^ + (ayj'] . (4) = ^ f^e'^^' /^^^ K,{ax' /Arl) 



(5) 



(6) 



(7) 



(8) 



Here, m is the effective electron mass, Swq is the char- 
acteristic confinement energy, 2d measures the center to 
center distance between the quantum dots, and a de- 
notes the ratio of the confinement strengths in the x and 
y directions. Moreover, we introduce the characteristic 
oscillator length rg = yJh/mujQ. The potential is shown 
in Fig. [1] together with a numerically calculated charge 
density. In the limit li ^ 0, the potential reduces to 
that of a single quantum dot. In our calculations we 
use material parameters typical of GaAs (m = 0.067me, 
£r = 12.9). We consider both the quasi-one dimensional 
limit a 3> 1 and the two-dimensional case a = 1. 

The exchange coupling between the two electrons is 
a purely orbital effect which arises as a consequence of 
the Pauli principle and the Coulomb interaction which 
lead to a splitting J = Ea — Eg of the lowest eigenvalue 
Eg corresponding to a symmetric orbital wavefunction of 
the two electrons, 4's(ri,r2) = 4'5(r2,ri), and the low- 
est eigenvalue Ea corresponding to an anti-symmetric 
orbital wavefunction, \E'yi(ri,r2) = — ^'A(r2, fi). Due 
to the Pauli principle the orbital part of a singlet state 
must be symmetric, while the orbital part of a triplet 
state must be anti-symmetric. The splitting of the orbital 
wavcfunctions may thereby be mapped onto an effective 



as the (regularized) Coulomb interaction in one dimen- 
sion. Here, Kq is the zeroth-ordcr modified Bcssel func- 
tion of the second kind. The exchange coupling can 
now be calculated using finite-elements by mapping the 
one-dimensional two-particle problem onto an effective 
two-dimensional single-particle problem: We consider the 
two-particle wavefunction '^(x\^X2) as describing a sin- 
gle fictitious particle with spatial coordinates r = (xi, xq) 
and momentum p = (pxj^,Px2)- Mathematically, the cor- 
responding singlc-particlc-likc Hamiltonian then reads 



H = p^/2m + W{r), 



(9) 



where 



Wir) = Vixi) + Vix2) + Ca^ilxi - X2I) (10) 



is the effective external potential that the fictitious par- 
ticle experiences. 

In this reformulation of the problem, the symmetry 
of the original two-particle wavefunction enters via the 
boundary condition along the diagonal X2 = xi. Sym- 
metric wavcfunctions fulfill ^'5(2:1,2:2) ~ ^'s(a;2,a;i) and 



consequently d^^^ s{xi, X2)\a 



dx2^sixi,X2) 
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(Neumann condition), while anti-symmetric wavefunc- 
tions fulfill '^a{xi,X2) ~ — ^^(a;2,xi) and thus 
1 X2)\x2=xi — (Dirichlet condition).— Since W^(r) 
is a confining potential, eigenfunctions go to zero in the 
limit |r| — > oo. In the numerical calculations we as- 
sume that the eigenfunctions are zero outside a certain 
finite range, and we check that the results converge with 
respect to an increase of this range. Thus, we only 
need to solve a one-particle problem on a finite-size two- 
dimensional domain with well-defined boundary condi- 
tions. This class of problems is computationally cheap 
with available finite-element method packagesJ^ 

Before discussing the numerical results we briefly 
review the standard approximations In the Heitlcr- 
London approximation the exchange splitting is cal- 
culated as Jhl = - with the 
Heitler-London wavefunctions |±) = (|L)i|_R)2 ± 
\R)i\L)2)/y/2{l ± |(L|i?)|2), where H is the full two- 
particle Hamiltonian, and \L) and \R) are the single- 
particle Fock-Darwin ground states of a single quan- 
tum dot centered at = (— d, 0) and r^j = (d, 0), 
respectively. The Heitler-London approximation can 
be improved by including doubly occupied spin singlet 
states and diagonalizing the Hamiltonian in the result- 
ing Hilbcrt space. This is known as the Hund-MuUiken 
approach and yields the expression Jhm = V — Ur/2 + 
^\/U^ + 16t2. Here, Ur and tr are the on-site Coulomb 
interaction and the tunnel coupling, respectively, renor- 
malized by the inter-dot Coulomb interaction as de- 
scribed in Ref. 0, while V (not to be confused with 
the confinement potential) is the difference in Coulomb 
energy between the singly occupied singlet and triplet 
states. Additional details about the approximative meth- 
ods are given in Appendix [XI 

If the inter-dot Coulomb interaction is negligible, the 
renormalizcd quantities [/,,, and t,. reduce to their bare 
values, U and t, while V = 0, and if moreover t/U ^ 1, 
the Hund-MuUikcn expression reduces to the standard 
Hubbard expression Jh = AP/U. The Hubbard approx- 
imation, which always predicts a positive exchange en- 
ergy, obviously cannot explain that the exchange energy 
with an applied magnetic field can be negative. This fail- 
ure can be corrected by retaining the inter-dot Coulomb 
interaction, and in the limit tr/Ur <^ 1, the Hund- 
MuUiken approximation then yields the extended Hub- 
bard approximation: Jj^ = At'^/Ur + V. The energy dif- 
ference V is important for the prediction of the exchange 
coupling at finite magnetic fields, where it allows for the 
predicted exchange coupling to become negative. 

In Fig. [2] wc show numerical results for the exchange 
coupling as a function of the inter-dot distance with 
different values of the confinement energy hujQ for the 
quasi-one dimensional case a = 10 3> 1. Together with 
the numerical results we show the Heitler-London, the 
Hund-MuUiken and different variations of the Hubbard 
approximations. The validity of the Heitler-London ap- 
proximation is strongly dependent on dimensionality due 
to the increasingly dominating Coulomb interaction in 




d/ro 

FIG. 2: (Color online). Exchange coupling as function of 
interdot distance in a quasi-one dimensional double quantum 
dot, a = 10. The distance d is measured in units of ro = 
yfi/rriZjo- Together with the numerical results we show the 
exchange couplings obtained with the Heitler-London Jhl, 
the Hund-MuUiken Jhm, the extended Hubbard J^, and the 
standard Hubbard Jh approximations. We also show J^ — V, 
where V is the Coulomb energy difference between the singly 
occupied singlet and triplet states. 

lower-dimensional systemSfi^ and for the quasi-one di- 
mensional case Jhl is negative in the entire range con- 
sidered for HajQ < 8 meV. The standard Hubbard ap- 
proximation predicts reasonably well the d-dependence, 
while both the Hund-MuUikcn and extended Hubbard 
approaches lead to (unphysical) negative values of the 
exchange coupling for a wide range of system parame- 
ters. We discuss these discrepancies in more detail when 
we consider the two-dimensional case below. Confine- 
ment energies larger than 18 meV are required for these 
approximations to yield positive exchange couplings for 
all interdot distances. For higher values of a, correspond- 
ing to stronger confinement in the y-direction, the range 
of validity of these approximations is further reduced. 

B. Two-dimensional case 

We next consider the two-dimensional case a = 1. In 
two dimensions the exchange coupling is strongly depen- 
dent on applied magnetic fields, and we include a mag- 
netic field perpendicular to the motion of the electrons 
by the substitution p p-l-eA, where A = Bz{—y,x)/2 
is a vector potential corresponding to the applied mag- 
netic field BzZ. The Zeeman term does not affect the 
exchange coupling and is trivial to include in final total 
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FIG. 3: (Color online). Exchange coupling as function of 
interdot distance in a two-dimensional double dot, a — 1. See 
Fig. [5] for details. The vertical lines denote the two values of 
d/ro for which the dependence on the magnetic field is shown 
in Fig. m 

energy calculations. 

Rather than mapping the two-dimensional two-particle 
problem onto an effective four-dimensional one-particle 
problem, we construct a two-particle basis from single- 
particle eigenstates (/>i(r) with eigenenergies e,; found 
by diagonalizing the single-particle Hamiltonian h{r) = 

''^^In '' — ^ ^i''^)- again using finite-element methods.~ 
The (unsymmmetrizcd) two-particle basis functions are 
then ^^ij (ri, = 0i(ri)0j (r2), in terms of which the 
matrix elements of the two-particle Hamiltonian read 

= (e, +£j)'5»,."5jj' + (*»j|C|*»'j')- (11) 

The Coulomb matrix elements are evaluated by inserting 
a set of two- particle states constructed from orthonormal- 
ized Gaussian single-particle wavef unctions. From the 
low-energy spectrum of H we then obtain the exchange 
coupling J. The details of this procedure are described 
in Appendix [Bl 

In Fig. [3] we show the results for the two-dimensional 
case a = 1. While the standard Hubbard approximation 
predicts well the d-dependence of the exchange coupling, 
the Heitler-London and the Hund-MuUiken approxima- 
tions yield predictions that in certain parameter ranges 
deviate significantly from the numerical results. In par- 
ticular, in the case Hloq = 4 meV a range of distances 
exists around d = rp , where both approximations predict 
negative exchange couplings. It is well-known that the 
Hcitlcr-London approximation fails at short distances. 



FIG. 4: (Color online). Exchange coupling as function of 
magnetic field in a two-dimensional double dot. Results were 
obtained with fkjo = 6 meV. See Fig. [2]for details. 

when the overlap of the Heitler-London wavefunctions 
becomes large, and that the range of validity is reduced 
as the ratio between the Coulomb and confinement en- 
ergy is increasedji^ This explains why the discrepancies 
are less pronounced in the case Hloq = 6 meV. We conjec- 
ture that the poor predictions by the Hund~Mulliken and 
the extended Hubbard approximations are mainly due to 
the Coulomb energy difference between the singly occu- 
pied singlet and triplet states, denoted V, overestimating 
the effects of the inter-dot Coulomb interaction at short 
distances {d ^ vq), leading to a too low (or even nega- 
tive) exchange energy. For large distances {d ^ 2ro), this 
overestimation decreases and a better agreement with the 
full numerics is obtained. In the figure we also show 
Jj^ — y which predicts well the exchange coupling, indi- 
cating that the effects of the inter-dot Coulomb interac- 
tion indeed seem to be overestimated. With larger con- 
finement energies this overestimation becomes less signif- 
icant, and a better agreement with the numerically exact 
results is found. 

In Fig. |4] we show numerical results for the exchange 
coupling as function of the magnetic field B with differ- 
ent inter-dot distances d. Together with the numerical 
results we again show the Heitler-London, the Hund- 
MuUiken and different variations of the Hubbard approx- 
imations. The results show that none of the approxima- 
tions predict well the dependence of the exchange cou- 
pling over the full range of magnetic fields for short dis- 
tances d < tq. For the Hund-MuUiken and the extended 
Hubbard approximations we again attribute the discrep- 
ancy to an overestimation of the effects of the inter-dot 
Coulomb interaction. For large distances this overesti- 
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mation is less pronounced, and a good prediction of the 
qualitative features is obtained. 

IV. CONCLUSIONS 

We have presented numerically exact finite-element 
calculations of the exchange coupling between electron 
spins confined in low-dimensional nanostructures. We 
have tested a number of approximations often encoun- 
tered in the literature by applying them to a simple dou- 
ble dot potential and found that they only predict well 
the exchange coupling in restricted parameter regimes, 
when compared to numerical exact results. While the 
approximative schemes may yield some insight into the 
qualitative features of the exchange coupling, we find it 
unlikely that they would suffice in the exchange coupling 
calculations for actual experimental structures and ex- 
periments, having seen how they may fail even in the 
case of a simple model potential. 

We thank A. Braggio, K. Flensberg, L. H. Olsen, A. 
S. S0rensen, and J. M. Taylor for fruitful discussions. In 
particular, we thank A. Harju for helpful advice during 
the development of our numerical routines. APJ is grate- 
ful to the FiDiPro program of the Finnish Academy for 
support during the final stages of this work. 

APPENDIX A: APPROXIMATIVE METHODS 



where b is the magnetic compression factor b = w/luq, 
If) is the zeroth-order Bessel function, erf (x) is the error 
function, and we have introduced the dimensionless dis- 
tance d —^ d/rQ. The prefactor c is the ratio between the 

Coulomb and confining energy, c, = -. — -w^- 

In the Hund-MuUiken approximation the ex- 
change coupling is calculated by diagonaliz- 
ing the two-electron Hamiltonian in the space 
spanned by *^^(ri,r2) = <I>±d(ri)$±d(r2) and 
*|(ri,r2) = [$+rf(ri)$_d(r2) ± $-d(rl)$+d(r2)]/^/2, 
where ^±d are the orthonormalized single-particle 

states $±d = {^±d - g(p^d)/\/^ - 2Sg + g'^, with 

g = (I — Vl — S^) I S. This leads to the expression 
Jhm = V- Ur/2 + + 16*2, where^ 

U = t - w ^ {^±d\h\<i>^d) - {"^llC \^^^) /V2, 
V^V--V+^{'f^C\'f^_)- (*^| C ) , 
Ur^U-V++X = C 1*^,,) - (*f I C ) 

+ {nd\C\^^d)- (A3) 

The Coulomb matrix elements are given by Burkard et al. 
in Ref. [2 and are applicable to any model potential for 
which the corresponding single-dot potential is a simple 
harmonic oscillator. Thus only the matrix element t is 
different for our model potential. We find 



In the quasi-one dimensional limit a 3> 1 we have 
evaluated the approximative methods numerically using 
Mathematica, ensuring convergence of the results with 
respect to a screening length S ^ of the regularized 
Coulomb interaction. In the following, we list analytical 
expressions obtained for the various approximative meth- 
ods presented in the article for the two-dimensional case 
a=l. 

The single-dot potentials corresponding to the double 
dot potential in Eq. (|4]) are those of a harmonic oscillator 
centered at (id, 0). The single-dot orbitals are thus the 
Fock-Darwin states shifted to (id, 0). For d ~ Q the 
Fock-Darwin ground state is 



(Al) 



where oj = ijo^/I + ujJJuj^ with ujl denoting the Lar- 
mor frequency ujl ~ eB/2mc. In the presence of 
a magnetic field given by the vector potential A — 
Bz{~y,x)/2, shifting the ground state to (±d, 0) adds a 
phase factorof e^*^'*/2'i3^ where Ib is the magnetic length 
Ib = \JhcleB. We thus obtain the single-dot orbitals 

V±d{x,V) = e^^y'^^'^^Bip(^x T d,y), where f±d{x,y) then 
denotes the single-dot orbital centered at (±d, 0). 

Using these single-dot orbitals we obtain for the ex- 
change coupling in the Heitler-London approximation 



sinh [2^2(26-1/6)] 
2d 



-bd- T 

e M 



{l-e-"^^} 



2d^ {l 



erf 



o(M2)} 
(Vbd)] 



(A2) 



hujQ 1 — 5*2 



^ (l - e-"'^') + rf2crfc (dVb) 



(A4) 



where erfc(a;) is the complimentary error function. 



APPENDIX B: NUMERICAL METHODS 



Here, we discuss the numerical method used in the 
two-dimensional case a = 1. We use finite-element 
methods to solve the single-electron problem given by 
the single-particle Hamiltonian h in Eq. (IT])J^ The full 
two-electron problem is then solved by expressing the 
two-electron Hamiltonian in Eq. ([T]) in a basis of prod- 
uct states of single-electron solutions \ipn), in terms of 
which the matrix elements are given by Eq. pT|) . To 
evaluate the Coulomb elements the single-electron eigen- 
states are expanded in an orthonormalized basis of 2D 
Gaussians 4'n^,ny{x, y) = x'^^y"'^e~^ 1"^ , where and Uy 
are positive integers or zero. The Coulomb matrix ele- 
ments between product states of 2D Gaussians can be 
determined analytically, and we state the result here for 



convenienc( 



.23 
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2 / , X «/2 Lni/2J L"4/2J 

S4=0 

r(7ii + i) r(n4 + i) 



r (si + 1) r (711 - 2si + 1) r (54 + 1) r {m - 254 + 1) 

r [{ni + ?i3 ~ 2si - 2s3 + 1) /2] r [(?i2 + 714 ~ 2s2 - 2s4 + 1) /2] 
r [(n - 2s) /2 + 1] 



X 2("-2''+i)/2pj(^_2s + l)/2], 



for rii + na and 712 + 714 even and zero otherwise. Here 

Cijkl = ('/'n,,.,n„,.0ri,,,,ri^,, |C|0„^_,,„^_,(/)„^_,,„^ ,) while 

r(x) is the Gamma function and [n/2j indicates floor- 
ing of half-integers. Above, we have introduced 77,1 = 
nx,i + rix.k, n2 = Uy^i + Uy^k, "-3 = J^rj + n^^i and 
77,4 = "yj + ny^i , while 77, = a-nd s = Z^i'**- 

The two-particle Hamiltonian matrix resulting from this 
procedure may then be diagonalized in the subspaces 
spanned by the symmetric and antisymmetric product 
states, respectively, to yield the exchange coupling. Be- 
cause the expansion in 2D Gaussians becomes increas- 
ingly inaccurate as the intcrdot distance is increased, we 
arc limited to intcrdot distances of the order of the char- 
acteristic oscillator length tq. The accuracy of the 2D 
Gaussian expansion at larger interdot distances could 
be greatly improved by using an expansion in relative 
coordinates 

The finite-element calculations of the single-particle 
states can be carried out with very high efficiency using 
existing finite-element packages^ and are not a limiting 
factor in terms of computational time or convergence. 
Also, the Coulomb matrix elements Cijki may be pre- 
calculated and saved in a lookup table, such that the 
largest portion of the computational time is spent as- 
sembling the two-electron Hamiltonian matrix. For each 
matrix element, a total of N"^ lookups in the Cijki table 
are required, where N is the number of 2D Gaussians 



included in the expansion set. A significant reduction 
in computational time is accomplished by utilizing the 
symmetry of the Hamiltonian in the product state basis, 
limiting the calculation to matrix elements which differ 
by more than a simple complex conjugation. For the re- 
sults presented in this article, a total of f 00 2D Gaussians 
were used to ensure that the results obtained may essen- 
tially be considered exact. With this basis set and a total 
of 7^ = 49 single-particle product states, the calculation 
of the exchange coupling takes approximately two and 
a half hours on a standard computer equipped with an 
Intel Core2 Duo 1.86MHz CPU. As few as 25 Gaussians 
arc in many cases sufficient to produce results that are 
within 10% of the exact results, and in that case a single 
calculation only takes about 5 minutes. 

The use of finite-element methods for solving the 
single-electron problem makes it easy to construct the 
two-electron Hamiltonian, even if analytic expressions for 
the matrix elements of the single-electron Hamiltonian in 
the basis of 2D Gaussians cannot easily be obtained. This 
makes the method very flexible, and only little work is 
required to solve problems with different choices of po- 
tentials. We have verified our numerical implementation 
against the results in Ref. d as well as for the simple prob- 
lem of two opposite spin particles in a two-dimensional 
parabolic potential, which can be solved analytically. 
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